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Abstract 

We study the efficiency of a two-channel reaction between two walkers on a finite ID periodic 
lattice. The walkers perform a combination of synchronous and asynchronous jumps on the lattice 
and react instantaneously when they meet at the same site (1st channel) or upon position exchange 
(2nd channel). We develop a method based on a conditional first-passage problem to obtain exact 
results for the mean number of time steps needed for the reaction to take place as well as for higher 
order moments. Previous results obtained in the framework of a difference equation approach are 
fully confirmed, including the existence of a parity effect. For even lattices the maximum efficiency 
corresponds to a mixture of synchronous events and a small amount of asynchronous events, while 
for odd lattices the reaction time is minimized by a purely synchronous process. We provide an 
intuitive explanation for this behavior. In addition, we give explicit expressions for the variance 
of the reaction time. The latter displays a similar even-odd behavior, suggesting that the parity 
effect extends to higher order moments. 
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1. INTRODUCTION 



In experimental systems there is always a finite discretization of the time window used for 
measurements. Thus, the degree of synchronicity for the dynamics of the different system 
constituents may vary when the size of the time window is changed, either intentionally 
or accidentally as a result of random fluctuations in the clock mechanism. But even for a 
fixed time resolution there may be intrinsic features of the experimental system that lead 
to the coexistence of synchronous and asynchronous events in the course of the evolution. 
Consider e.g. diffusion of an ensemble of identical particles in a disordered or a randomly 
fluctuating medium 0, 0| • The diffusivity of a particle will then be different depending on 
its location, and this will lead to the coexistence of different characteristic time scales in the 
system, resulting in partial desynchronization for the dynamics of its constituents. 

Recently, the author and coworkers used a simplified lattice version of this problem to 
study the role of synchronicity on the efficiency of a diffusion-controlled two-channel reac- 
tion between two particles 0]. Here, the parameter chosen for the characterization of the 
reaction efficiency is the mean number of time steps to reaction. The effect of synchronic- 
ity was investigated by studying how this quantity is affected by the interplay between the 
transport dynamics and the geometric characteristics of the lattice, i.e. size, dimensionality 
and boundary conditions. Finite size effects were found to play an important role, thus 
emphasizing the relevance of geometric constraints for systems with a small number of con- 
stituents. On the other hand, a variety of interesting problems in statistical physics may be 
recast in terms of a random walk in a lattice of a small size. An example thereof is a family 
of ruin problems where the capital of each gambler is typically a small number 0, 0|. 

In ref. Q the analytic results for the ID case were based on a difference equation 
approach, whereas numerical simulations were used in higher dimensions. However, it is 
also instructive to compute the mean reaction time for this problem by other methods, 
e.g. approaches based on the theory of finite Markov processes (HI 0, 0| and on generating 
function techniques 0. The first method has the advantage of being applicable in higher 
dimensions and to systems lacking translational invariance (e.g. with reflecting boundaries), 
while the second emphasizes the correspondence between diffusion-controlled reactions and 
first-passage problems In this paper, we shall apply the latter to determine the mean 
number of time steps necessary for the reaction to take place as well as the variance of this 
quantity. An important advantage of this approach is that the variance and higher order 
moments can be computed by a straightforward differentiation of the relevant generating 
function, as opposed to the method of difference equations, where moments are coupled to 
each other via a hierarchy of equations |3| . 

The work plan is as follows: in Section|2]we introduce the model and its reduced represen- 
tation in a comoving reference frame. Section El introduces the generating function approach 
to compute the reaction time and higher order moments. For completeness, we first present 
the standard formalism to deal with some simple cases where only a single reaction channel 
is active. We then treat the general case, which is solved in two stages. First, we show that 
our model is equivalent to a model with a relaxed collision rule and then we solve for the 
latter. Section H] contains the main results for the reaction time and its variance. Finally, 
Section El summarizes the main conclusions. 
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2. THE MODEL 



The system to be studied consists of two coreactants A and B performing symmetric 
nearest neighbor random jumps on a iV-site periodic lattice at discrete time steps. The 
reactants, also termed walkers in what follows, are assumed to react with each other whenever 
they meet at the same lattice site or attempt to exchange positions. Each of such "collisions" 
results in an instantaneous irreversible reaction. Regardless of the particular outcome of the 
reaction , the mean reaction time can be identified with the mean number of time steps 
elapsed until the collision takes place; the larger this quantity, the less efficient the reaction 
will be. Each time step will be considered to be an elapsed time unit regardless of the lattice 
size. In the course of the dynamics, the following joint events may occur: 

1. with probability p both walkers hop simultaneously to randomly chosen nearest neigh- 
bor sites (synchronous event). 

2. with probability 1—p one of the walkers (no matter whether A or B) hops to a nearest 
neighbor site while the other one remains immobile (asynchronous event). 

Thus, the characteristic parameter p interpolates between the asynchronous case (p = 0) 
and the case of two simultaneously moving walkers (p = 1). Occasionally, we shall refer 
to the limiting cases p = and p = 1 as "the purely asynchronous case" and "the purely 
synchronous case" respectively. The purely asynchronous case was studied by Montroll 
G3 and Montroll and Weiss in one, two and three dimensions. For the particular 
case of a periodic ID lattice they obtained: 



where (n) is the initial-condition-averaged mean reaction time. Note that in this case the 
only active reaction channel is same site occupation (SSO), since reaction by position ex- 
change of the walkers, i.e. by nearest neighbor crossing (NNC) is not possible. 

However, in the purely synchronous case, both reaction channels SSO and NNC will be 
active if the total number of sites N is odd, otherwise only one channel will be available, 
depending on the initial location of the walkers. A consequence is that the analytical ex- 
pression for (n) in terms of the size of the ID lattice depends on the parity of N 0, . 
This even-odd effect translates mathematically as follows 



In contrast to the above cases, as soon as < p < 1, reaction by NNC becomes possible 
regardless of the value of N and the initial two-walker configuration. 

Due to the translational invariance of the lattice, only the relative motion of both walkers 
is relevant for the computation of the collision time. This has two important consequences. 
First, it tell us that the physical distinguishability of the walkers is irrelevant for the solution 
of the problem. To emphasize this, the walker labels A and B have been left out in the scheme 
displayed in fig. The second consequence is that it is convenient to choose the reference 



1 The case A + B — > inert, is of great historical interest, see e.g. refs. 1 1 (I 111! ll'A \13L li t ll 




N{N + l)(N + 2)/(12(N-l)) 
(JV + l)(iV + 3)/12 



for N even, 
for N odd. 



(2) 
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frame in such a way that one of the walkers (say B) is at rest. In this comoving frame, 
walker A will either hop to a nearest neighbor (with probability 1 — p), hop to a next to 
nearest neighbor site (with probability p/2) or remain immobile (again with probability p/2) 
as a result of the rules 1 and 2 prescribed above. In this one-walker representation, walker 
B plays the role of a stationary trap T, as indicated in fig. ^p; any time walker A reaches 
or attempts to overcome the site at which walker B is placed, the instantaneous reaction is 
triggered and the dynamics is stopped. In what follows, we shall therefore refer to walker B 
as "the trap" and to walker A as "the walker" when we work in the comoving frame. 

For convenience, let us place the origin (site 0) of the comoving frame at the initial 
position of walker A and then number the remaining sites, say clockwise, from 1 to N — 1. 
Denoting by jr the coordinate of the site at which the trap is located, the distance in lattice 
spacings between both walkers will be d = min(j<r, N — jx)- 

3. GENERATING FUNCTION APPROACH 

Our next goal will be to derive an expression for the mean number of time steps (n)j T 
necessary for the reaction to take place for a given value of the coordinate jx characterizing 
the initial condition. As a starting point, we take the equations that would govern the 
sojourn probabilities of an unrestricted walk if there were no interaction between the walker 
and the trap. These equations read 

V n+ i(j) = i [Pn(j - 2) + 2V n (j) + V n (j + 2)] + l^l [V n (j - 1) + V n (j + 1)] (3) 

where V n (j) is the probability to find the walker at a given site j after n time steps (j may 
take integer values from and N — 1 and the site addition and subtraction is performed 
modulo N). The first term in the r.h.s. of eqs. Q is the contribution due to the synchronous 
events, by which the walker either remains at rest or it moves two lattice sites clockwise 
or anticlockwise. The second term describes jumps by one lattice site yielded by the asyn- 
chronous events. In accordance with our definition for the origin, eqs. must be solved 
using the deterministic initial condition P = (V (0), V (l), • • • , V (N — 1)) T = (1,0,..., 0) T . 

3.1. A simple case 

The next step is to incorporate the walker-trap interaction to the above formalism. Let 
us first consider the situation where reaction by NNC is precluded. This holds if and only 
if TV and the walker-trap separation d are even integers. In this case site jx can be viewed 
as a reactive site, also termed "r-site" in what follows. Clearly, the mean reaction time will 
be given by the mean first-passage time {n) - t of the walker at site jx- 

oo 

K T = $>.F n (j T ), (4) 

n=l 

where ^F n {j) is the probability of visiting a given site j for the first time after n time steps. 
Eq. (HJ) can be expressed as follows 

(5) 

2 = 1 



d 
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where T{j,z) = Y^n=o-^n{j) z n is the generating function of the first-passage probabilities 
T n {j) and the limit 2; — 1 is taken from below. On the other hand, this function can be 
directly related to the generating function V(j, z) = ^2^=0 V n (j) z n . To do so, one uses the 
fact that the sets of probabilities V n (j) and T n {j) are linked to each other via the equation 



(6) 



k=l 



The discrete convolution in the r.h.s. of eq. (jUJ) is equivalent to the product V(j, z) T{j, z) 
in the reciprocal generating function space. Thus, we have 



n=0 



V{0,z) 



The collision time (n)j T then follows from eqs. © and ((7j): 



(n) 



d V(jt,z) 
jT dz V(0,z) 



(7) 



(8) 



2=1 



3.2. General case 

Let us now extend these results to the general case where N, d and the synchronicity 
parameter p take arbitrary values. The strategy to tackle the problem will be as follows: 
we shall not deal with NNC events directly, but rather introduce a model with a relaxed 
definition of collision, show its equivalence to the original one by virtue of the topological 
restrictions imposed by the ID lattice and then compute the reaction time for the relaxed 
model. For brevity, let us respectively refer to the original and the relaxed models as "model 
I" and "model II" in what follows. 

In model II, one assumes that the walkers react instantaneously by SSO or when they 
jump to nearest neighbor sites, i.e. by nearest neighbor occupation (NNO). The represen- 
tation of the system in the comoving frame will consist of a walker in a lattice with three 
r-sites, as shown in fig. |21 for a 7-site system. According to the site numbering introduced 
in Sec. 121 these three r-sites will have the coordinates jx (where the immobile reactant is 
located), ]t — 1 and jt + 1- By definition, the walk will automatically terminate when the 
walker lands on any of the three r-sites. For the purpose of computing (n)j T , it is easy to 
realize that the system can be unfolded into an equivalent system with a non-periodic lattice 
by introducing an additional fictitious r-site, as shown in fig. El for N = 7. 

On the other hand, let us again consider the walker-trap representation of the model I 
depicted in fig. [TJd for iV = 7. Along the same lines as above, the system can be unfolded 
into an equivalent transformed lattice with two trapping sites T and N — 1 = 6 non-trapping 
sites (cf. fig. HJ). Each site T can then be replaced with two fictitious r-sites, as shown in fig. 
13] 2 . Thus, we realize that model I embedded in a periodic lattice with iV sites is equivalent 
to model II in a periodic lattice with iV + 2 sites. 



2 If the dynamics is purely asynchronous (p = 0) , only one r-site at each end will be needed, since jumps 
by two sites are not possible in this case. 
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Next, let us compute the reaction time {n)j T for model II. The situation is now more 
complex than the one described in the previous subsection, since we do not know a priori 
at which of the three r-sites Jt — 1, Jt or jt + 1 the reaction will occur. Yet, it can still 
be formulated as a (conditional) first-passage problem. The key quantity in this case is 
T n {s | s x ,n x ; s 2 ,n 2 ), i.e. the probability that a random walker arrives at site s for the first 
time at the n-th time step after having visited the sites s x and s 2 exactly ri\ and n<i times 
respectively. For a given initial condition with a fixed value of jx, the mean collision time 
is now 



Ht = ^{^n(jT|jT-l,0;jr + l,0) 

n=0 

+M3T ~ 1 1 Jt, 0; j T + 1,0) + T n {3T + 1 1 j T - 1, 0; JT, 0)} . 
The r.h.s. of eq. can again be expressed in terms of the generating functions 



(9) 



(n) 







{^(Jt|jt-1,0;jt + 1,0;z) 



/; dz 

+F{j T - 1 I Jt, 0; Jt + 1, 0; z) + F(j T + 1 1 j T - 1, 0; j T , 0; z)} 



where 



z=l 



(10) 



F{s | s x ,n x ; s 2 ,n 2 ; z) = y^njs \ s x ,n x ; s 2 ,n 2 ) z n . 

n=0 

To compute this generating function, let us first observe that 



(11) 



V n (s | si,0; s 2) 0) = y^^fe(.s | si,0; s 2 , 0)P n _ fc (0 | si - s, 0; s 2 - s,0), s^0, si, s 2 , (12) 



fc=i 



where P n (s | s x , 0; s 2 , 0) is the probability that the walker is at site s at time n conditioned 
to its not having visited sites s x and s 2 but regardless of any previous visits to s. This 
probability must fulfil the initial condition Vq(0 \ s x — s, 0; s 2 — s, 0) = 1. A simple calculation 
shows that eq. (|T2*|) leads to the relation 



T{s | s x ,0;s 2 ,0;z) 



V{s | s 1 ,0;s 2 ,0;z) 
V(0 | s x — s, 0; s 2 — s, 0; z) 



(13) 



where 



V{s | s x , n x ; s 2 , n 2 ; z) = ^ I s i> n i5 s 2, na) ^ 



(14) 



n=0 



This generating function can be computed in terms of V(s, z). For details we refer to the 
book by Weiss (ref. [H, chapter 4). The result for n x = n 2 = is 



2 ^ 2^ fe (o) 

V(s | s 1; 0; s 2 , 0; *) = V(s, z) - ^ -^y- P(s - s fc , z), s ^ si, s 2 



(15) 
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where T>k(0) and T>(0) are the following determinants: 



P(0) = 

P!(0) 

V 2 (0) 



V{0,z) V{s 2 -s u z) 
V( Sl -s 2 ,z) V(0,z) 

V(s 1 ,z) V{s 2 - sx,z) 
V(s 2 ,z) V(0,z) 

V{0,z) V{s x ,z) 
V(s 1 - s 2 ,z) V(s 2 ,z) 



(16) 
(17) 
(18) 



We can now use these expressions to compute explicitly the mean collision time {n)j T - From 
eqs. (fTUJl and (JT3Jl we have 



(n) = d f V(j T \j T -l,0;j T + l,0;z) V(j T - 1 | j T , 0; j T + 1, 0; z) 
WjT dz \ V(0 1 1, 0; iV — 1, 0; z) V(0 | 1, 0; 2, 0; z) 

, V(j T + l\j T -l,0;j T ,0;z) 



V(0\N- 1,0; N- 2,0; z) 



(19) 



2 = 1 



(recall that the arguments of the different generating functions are evaluated modulo N). 
Thus, the reaction time can again be expressed in terms of V(j,z). This quantity is easily 
computed from the relation 



N-l 



(20) 



j'=0 

where the ]3(j)'s are the single step probabilities for the random walk. Next, let us introduce 
the Fourier transform of the single step probabilities: 



P 



2tt/ 



N-l 



P(k) exp I 

fc=0 ^ 



/ 2nik I 



N 



(21) 



This expression can be used to represent V(j, z) as follows 



N-l 



k=0 



exp (2irikj/N) 
zp(2nk/N) ' 



P 



(22) 



(23) 



where p(-) is computed from the following single step probabilities: 

p(2)=p(-2) = |, p(l)=p(-l) = ^, p(0) 2 . 

We now have the necessary ingredients to evaluate the functions V{j\ \j 2 ,0;j 3 ,0; z) ex- 
plicitly and to compute {n)- T and (n) via the relations (fT3j) and (fl9*)l. The initial-condition- 
averaged reaction time (n) is obtained by averaging over a uniform ensemble of all possible 
non-reactive configurations characterized by distinct values of jt, i.e. 



N-2 



(n) 



N-3 



(24) 



JT=2 
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3.3. Higher order moments 



An important advantage of our generating function approach is that once T(i \ j, 0; k, 0; z) 
is known, the computation of higher order moments can be carried out straightforwardly by 
deriving with respect to z. Indeed, one has 



(n m >. 



JT 



z^) {F(j T \j T -l,0;j T + l,0;z) 



+HJt - 1 I Jt, 0; J T + 1, 0; z) + T{j T + 1 1 j T - 1, 0; j T , 0; z)} 
In particular, the second-order moment 



2=1 



(25) 



'%t = J2 n2 {^Ut I Jt - 1, 0; j T + 1, 0) + ^.(j't - 1 | j T , 0; j r + 1, 0) 



n=0 
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+^ n (j T + 1 1 j T - 1, 0; j T , 0)} = ( z — ) {^0t I Jt - 1, 0; j T + 1, 0; *) 



the variance 



+F{j T - 1 1 j T , 0; j T + 1, 0; z) + T{j T + 1 1 Jt - 1, 0; j T , 0; z) } 



and its average over the different initial conditions 



N-2 



are easily obtained. 



RESULTS 



4.1. Behavior of the reaction time 



2 = 1 



(26) 



(27) 



(28) 



In order to obtain the results for model I, the lattice size N must now be decreased by two 
units in the expressions for (n) and (v) obtained in the framework of the above formalism. 
The expressions for (n) as a function of p for small values of iV 3 are displayed in Table UJ 
It is seen that these are ratios of two polynomials whose complexity grows with increasing 
lattice size. These results are in full agreement with previous findings for (n) obtained in 
the framework of a difference equation approach by the author and coworkers Fig. El 
displays the (n)- plots computed from the ratios of polynomials given in Table E for small 



3 For large values of iV, the evaluation of the generating function becomes rather time-consuming, except 
in the limiting cases p — and p = 1. 
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values of the lattice size N. Let us now summarize their main properties and provide an 
intuitive explanation for the observed behavior. 

The plots confirm the validity of the expressions (0) and (J2J) for the two limiting cases 
p = and p = 1. For a fixed value of p the collision time always increases with the lattice 
size. On the other hand, the behavior of (n) as a function of p is not always monotonic. For 
N = 2 (the smallest physically interesting lattice) the most effective process is the purely 
asynchronous one, and (n) increases monotonically with p. For N = 3 the reaction time 
does not depend on p. For N = 4, 6, 8 the reaction time decreases with increasing p in a 
wide regime of p values, but it then increases again for a sufficiently large value of p. In 
contrast, for N = 5,7 and 9 (n) decreases monotonically with p. 

These results as well as numerical simulations for larger lattices confirm the existence of 
an even-odd effect for N > 3. For odd values of N, (n) decreases monotonically with p, 
whereas for even values of N it decreases up to a value p m in and it increases monotonically 
beyond this value. The value p min is rapidly shifted to one as N increases, and the curves 
get closer and closer to the linear law (n) = a{l + p) (with a > 0) predicted by the 
continuum approximation 3]. Deviations from this behavior for small lattices, where the 
purely synchronous process is not always the most efficient one, may thus be regarded as 
an indication of the important role played by finite-size effects as well as by the spatial 
discretization imposed by the lattice. 

The even-odd effect described above is reminiscent of the one observed in the limiting 
case p — 1 (f| [l^]. In the case p = 0, only one of the reaction channels (SSO) will be 
open regardless of the value of N. In contrast, when p = 1 one has two distinct behaviors 
depending on the parity of N. For odd values of N, both channels are open. However, 
for even values of N reactions occur via a single channel for a given initial condition, i.e. 
through SSO when the distance d is even or by NNC otherwise. Now, when p is no longer 
strictly equal to one but still close to this value, both reaction channels should be open, but 
the above parity effect still holds in a statistical sense, since NNC (SSO) reactions will be 
rare for an even (odd) value of d. Hence the parity-dependent behavior of the curves (n) (p) 
in this regime. 

The non-monotonic behavior of (n) for even values of N can be explained in terms of a 
competition between synchronous and asynchronous events. While synchronous transport 
is the most efficient mechanism to bring distant walkers in the vicinity of each other, it 
does not always maximize the efficiency of the reaction once the walkers are within the 
typical interaction radius. In a prereactive configuration with the walkers sitting at next to 
nearest neighbor sites, reaction within the next step can only take place via a synchronous 
event. However, if they are at nearest neighbor sites, reaction through an asynchronous 
event becomes possible and it is then more efficient than reaction through a synchronous 
event (the reaction takes place 2 times out of 4 vs. 1 time out of 4). Thus, if d and iV 
are such that prereactive configurations with contiguous particles are favored, the resulting 
competition between synchronous and asynchronous dynamics will lead to an antiresonance 
of the reaction time (maximal efficiency) for a value p m in between and 1. The statistical 
weight of such configurations is apparently stronger in the case of an even lattice, thus the 
antiresonance effect still prevails after averaging over the initial conditions, leading to the 
observed non-monotonic behavior of (n). In the limit of large N, the role of diffusional 
transport becomes increasingly important for the efficiency, thus one has p m in —> 1- On the 
other hand, one can also show that the parity effect is washed out in the diffusive limit, 
where the SSO and the NNC reaction channels become indistinguishable 0]. 
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4.2. Behavior of the variance 



The analytic expressions for (v) as a function of N in model I are given in table |H] The 
p-behavior of these functions follows essentially the same law as (n) (cf. figEJ). For odd N 
the behavior is monotonically decreasing, while for even N a minimum is observed which 
rapidly shiftes to p — 1 with increasing N. This suggests that the behavior observed for (n) 
extends to higher order moments; note, however, that the value p m in is slightly smaller than 
the one obtained from the (n)-curves, except for N = 4, where it turns out to be the same 
{Pmin = 2/3). Thus, (v) and (n) cannot be simultaneously minimized except in this case. 

In the limiting cases p = and p = 1 it is possible to simplify the expressions of the 
relevant generating function and thereby obtain general expressions for (v) for arbitrary 
lattice sizes. In the purely asynchronous case one has 

= N(N + im-2)(N + 2) 
w 30 V ; 

This expression is also recovered by taking the average of Montroll's original result for the 
variance over all possible initial walker-trap separations d [13]. On the other hand, the case 
p — 1 yields 

J N(N + 1)(N + 2){N 2 + 2N + 2)/(120(JV - 1)) for N even, 
W \ (N + l)(iV + 3)(iV 2 + 2N- 5)/120 for N odd. [ } 

The standard deviation a = ^/(v) thus turns out to be comparable to (n) in both cases. 
According to eqs. (J2HJ an d (JHOJ), one has respectively a oc N 2 /y/30 and a oc iV 2 /(2\/30) 
for large N, while eqs. ((TJ) and (J2J) give respectively (n) oc iV 2 /6 and (n) oc N 2 /12. The 
standard deviation and the mean value remain comparable for intermediate values of p. This 
is not surprising in view of the large variability characteristic of first-passage problems. 



5. CONCLUSIONS 

We have used a generating function approach to compute the mean reaction time be- 
tween two walkers performing a combination of synchronous and asynchronous jumps. The 
walkers react via two channels, i.e. same site occupation or position exchange. The reac- 
tion time and its variance display a different behavior for even and odd lattices, i.e. they 
behave monotonically as a function of the synchronicity parameter p for odd lattices and 
display antiresonances in the even-lattice case. This behavior has been explained in terms 
of a competition between synchronous and asynchronous dynamics. While the former favors 
diffusional transport over long distances, it may also lead to a decrease of the cross section 
once the particles are within the typical interaction radius, since they may avoid each other 
more easily when they hop simultaneously. The even-odd effect tends to vanish with in- 
creasing lattice size and in the continuum limit, since it is a signature of the discrete nature 
of the support. 

As pointed out in the Introduction, the interest of considering the effect of synchronicity 
on the behavior of a given system arises in a variety of different contexts. A fundamental 
motivation is provided by the fact that, while free particles move simultaneously under 
the action of a natural law at the microscopic level of description, in many real systems 
the existence of geometric and energetic constraints (e.g. activation energies) results in an 
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intrinsic nonzero degree of asynchrony in the dynamics of the constituents at mesoscopic 
time scales. Since in this context asynchronous dynamics can be understood as an expression 
of such constraints, it then becomes natural to ask how it affects the efficiency of a given 
physical process. In view of our results, the answer to this question cannot be considered to 
be straightforward. Moreover, it may turn out to be rather counterintuitive. 

Admittedly, the results depend strongly on the specific definition of the collision rules. 
Yet, the equivalence between models I and II suggests that certain features of the observed 
behavior might be characteristic of a class of small systems. Similar antiresonance phe- 
nomena in the encounter time have e.g. been recently observed in a model for target site 
localization of a protein on DNA [20]. On the other hand, the model shows that classical 
techniques inspired in first-passage problems can be successfully used to compute character- 
istic reaction times for complex processes involving more than one interaction channel. 

Possible extensions of the model include the higher dimensional case 4 , the case of more 
complex media and the evaluation of other quantities such as survival proba- 

bilities jiil . I25I ] , nearest-neighbor distance to the t rap |26| and concentration decays |27| in 
the framework of the many-particle problem [2^, |2£| or more ellaborate reaction schemes 
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N (n) 

2 2/(2 -p) 

3 2 

4 (10/3)(3p-4)/(p 2 + 2p-4) 

5 4(2p- 5)/(p 2 -4) 

6 (28/5) (p 2 - lOp + 10)/(p 3 - 4p 2 - 4p + 8) 

7 (4/3) (p 2 + 8p - 14) / (p 2 - 2) 

8 (12/7)(13p 3 + 6p 2 - \2<op + 112)/((p - 2)(p 3 + 6p 2 - 8)) 

9 10(2p 3 - 5p 2 - 16p + 24)/((p 2 + 2p - A)(p 2 - 2p - 4)) 

10 (22/9) (7p 4 - 7Qp 3 + 16p 2 + 288p - 240)/(p 5 - Qp 4 - I2p 3 + 32p 2 + 16p - 32) 

TABLE I: Analytic expressions for (n) for different lattice sizes in model I. 



N {v) 

2 2p/(2 - p) 2 

3 2 

4 (2/3)(192 - 316p + 152p 2 - 15p 3 )/(p 2 + 2p - 4) 2 

5 4(-2p 3 + 25p 2 - 80p + 84)/(p 2 - 4) 2 

6 (28/5)(-1096p + 768p 2 - 206p 3 + 28p 4 - p 5 + 512)/(8 - Ap + p 3 - 4p 2 ) 2 

7 (4/3)(-320p + 88p 2 + 8p 3 + p A + 252)/(p 2 - 2) 2 

8 (12/7) (-13p 7 + 252p 6 + 546p 5 - 4736p 4 - 2184p 3 + 35520p 2 - 50848p + 21504)/((p 3 + Qp 2 - 8) 2 {p - 2) 2 ) 

9 2(-10p 7 + 253p 6 - 816p 5 - 2748p 4 + 11872p 3 + 5552p 2 - 42496p + 29568)/((p 2 - 2p - 4) 2 (p 2 + 2p- 4) 2 ) 

TABLE II: Analytic expressions for (v) for lattice sizes up to N = 9 in model I. 
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a) 



b) 



FIG. 1: a) Two-walker system on a seven-site periodic lattice. Both walkers are represented by black 
circles. For convenience, the walker labels A and B have been left out (see text for explanation), 
b) Equivalent one-walker plus trap system. The trapping site is denoted by "T". The arrows in fig. 
QJa indicate that both walkers perform a synchronous step. In fig. QJ> this corresponds to a two-site 
jump of the walker. 
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FIG. 2: Representation of model II in the comoving frame. 
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FIG. 3: Lattice transformation for model II in the comoving frame. A periodic lattice with three 
reactive sites is equivalent to a non-periodic one with four r -sites 
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FIG. 4: Lattice transformation for the one-walker representation of model I displayed in fig. 
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